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I. INTRODUCTION 


Solutions to the equations of General Relativity that describe space-times containing 
closed timelike curves (CTC) have attracted significant interest as they revealed at least 
hypothetical theoretical possibility of travelling in time. Since the renowned publication by 
Kurt Gbdel fl] a number of causality violating solutions in GR as well as in modified theories 
of gravity have been constructed, among which we can name the Tipler-Van Stockum time 
machine generated by axially rotating distribution of particles 0. 0. the Morris-Thorne- 
Yurtsever transversablc wormhole M, the Gott time machine based on moving conical 
defects jB], the Ori dust solution [7], and the solutions in f(R) theories of gravity [8] and 
theories with non-minimal matter-curvature coupling [|9]. 

All questions about physics of time machines that could be posed in principle fall into 
three general categories: 

• Is there a physical way to create a time machine? 

• Is there any time machine solution that can be stable? 

• What dynamical behaviour would a physical system experience evolving in a time 
machine background? 

None of the questions have yet received a definite answer. 

The answer to the first question is believed to be negative. Extensive analysis of particular 
time machine solutions has demonstrated that in order to create a space-time with CTC 
one needs matter that violates strong, weak or null energy conditions of General Relativity 
(different solutions require violation of different energy conditions), and only eternal time 
machines can exist PEI EH- However we can not be sure that all matter in the Universe obeys 
these conditions. For instance, there are a number of models of the dark energy violating 
the null energy condition H2HI5], and this provides a way to by-pass the no-go statement. 

The second question was raised by Hawking in [15], where he conjectured that a space- 
time with CTC can be stable only on classical level, but will be unavoidably destroyed by 
quantum fluctuations of the metric. The real universal proof or refutal of the conjecture 
can be obtained only within a framework of a complete theory of quantum gravity. String 
theory opened a possibility to check the chronology protection condition in specific cases. 
In pTTj authors have shown that appearance of closed timelike curves in a certain (O-plane) 
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orbifold background would cause a Hagedorn transition that restructures the space-time 
transforming it into a chronologically safe configuration. So this result can be considered as 
a very accurate and nice supporting evidence in favour of the Hawking conjecture. On the 
other hand in jTHHSO] , it was demonstrated that the Godel type solutions can be smoothly 
embedded in the context of string theory. Closed timelike curves in that case are hidden 
behind the so called holographic screens and do not violate causality in the rest of the 
space-time. Thus the chronology is protected, but structure of the CTC remains unbroken 
by quantum effects. An intriguing observation has been made by authors of |2T] and [22], 
that from the point of view of the AdS/CFT correspondence, the existence of CTC in the 
bulk can be related to negative or exceeding one fermionic probability in the boundary held 
theory. 

The third category includes various types of the “grandfather paradox”. For a clas¬ 
sical wave equation on a non-globally hyperbolic space with CTC the possibility of self- 
consistent dynamics was demonstrated in [23] (2Tj. Classical mechanical billiards and their 
self-consistency conditions have been studied in [25] . A basis of states of a free quantum 
held theory in the Gott time machine has been constructed in [26], where it was shown that 
the causality violation leads to an emergence of an effective non-unitary interaction in the 
theory. Non-unitarity of interacting held theories in time machines was analyzed in (27]. 
Some authors even argued that evolution of a physical system along closed timelike curves 
can be studied experimentally by mean of simulation of emergent gravity in metamaterials 
[28], or a qubit interacting with an older version of itself [29]. 

However the question about properties of an interacting quantum held theory in a time 
machine background remains open, though the real “grandfather paradox” can take place 
only in a self-interacting system. When the notions of time ordering and unitarity are 
absent from the very beginning, it is unclear how to formulate an interacting held theory. In 
this paper we address this problem and by use of the AdS/CFT correspondence provide a 
constructive solution to it. Here we stand on the position that even if the presence of CTC 
causes breaking of unitarity in the boundary held theory [22], it should not be regarded as a 
big problem as long as we can formulate a prescription how to solve the theory. When one 
is trying to get an insight into physics of paradoxical systems, it is not very useful to rely 
on the “common sense” intuition and corresponding fundamental principles. 

The AdS/CFT correspondence provides an elegant way to address the paradox. In the 
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large N limit it relates quantum field theories to a classical gravity, and thus we can study 
properties of a quantum theory in the CTC background just by careful analysis of the dual 
Riemannian geometry, without any need to formulate special quantization rules that would 
be valid in the case of broken causality. 

The paper is organized as follows. In the next section we introduce a simple set up 


for the time machine in AdS 3 . In Sec. Ill we discuss geodesic structure of the spacetime, 
and suggest that it could lead to non-trivial effects in the boundary held theory. In Sec. 


IV we introduce a notion of timelike quasigeodesics that will be then used for connecting 


timelike separated boundary points. Finally, in Sec. [V]we formulate a precise algorithm for 
the Green function evaluation, provide the results of numerical simulations, and discuss the 
related phenomenology. 


II. TIME MACHINE IN AdS 3 

The eternal time machine solution in AdS has been suggested by Gott and DeDeo in [30] 
(for similar solutions containing CTC but collapsing into a BTZ black hole see [31]). Here 
we briefly recall its structure closely following the original text. 

The three dimensional global anti-de Sitter space-time can be thought of as a hypersurface 

- X „ 2 - A'J + X? + X% = -1, (1) 

embedded in a four-dimensional hat M 2,2 space-time with a metric: 

ds 2 = -dXl - dXl + dX\ + dX 2 . (2) 

In the Schwarzschild coordinates the embedding formulas are 

X 0 = \/l + R 2 cos t , (3) 

X 3 = vT + R 2 sin t , 

X 1 = R cos 0, 

X 2 = R sin 0 , 

where R G ( 0 , 00), t G (—00, 00), 0 G [ 0 , 27 t). 

The induced metric is then 

ds 2 = — (1 + R 2 )dt 2 ————r + i? 2 d0 2 . (4) 

1 + R z 
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FIG. 1 . Two equivalent unfoldings of the AdS^ spacetime with a conical defect a = The 

larger part of the spacetime outside of the faces is to be cut out, and the faces are identified. 
The only physical space is the narrow region between the faces. To construct a time machine the 
twisted unfolding is more convenient to use. 


A massive particle put into a three dimensional space-time removes a wedge with an angle 
deficit proportional to the mass of the particle, and edges (faces) emerging from this point¬ 
like particle. Points on the opposite edges of the wedge are identified, and the resulting 
space-time contains a conical defect, Fig0 a). When we are looking at the unfolding of 
the conical defect, coordinate locations of the edges do not have an independent physical 
meaning, and we are free to rotate them preserving the angular deficit. For our purposes it 
will be convenient to make the cut out “pizza slice” twist in time with a constant angular 
velocity in the reference frame of the massive particle, making a full rotation in a period 2n, 
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Figjl|b). Then for the trailing and leading faces of the wedge in the embedding coordinates 
we get: 


Xq = y/l + R 2 COS t x l 0 = yjl + R 2 COS t , 

X 3 = \/l + R 2 sin t X l 3 = y/l + R 2 sin t , 

X\ = Rcos(t — a/2) X[ = Rcos(t + a/2 ), 

X\ = Rsm(t — a/2) X l 2 — Rsin(t + a/2). 

Here a is the angular deficit of the conical spacetime. Integrating the spacetime stress-energy 
tensor over the angle, we can deduce that the effective mass concentrated in the interior of 
the bulk is 


M 


1 a 

8 G + 16v tG' 


(5) 


Here the second term is the mass of the point-like source, and the first one is the contribution 
from the negative AdS curvature. 

For a single static conical defect we can not make its angular deficit a larger than 2n. If 
the mass of the point-like particle exceeds the limit a = 2n (M > 0), the resulting space-time 
will be rather a BTZ black hole instead of a naked conical singularity |35J. 

However, if we boost the massive source, the removed wedge is effectively getting 
“squeezed” from the point of view of an external observer at rest (see Fig. [2](a)) . This 
provides a room for a second conical defect with a deficit angle (3 such that a + /3 > 2i r. 
In other words, relative motion can support the system of conical defects, preventing them 
from collapsing into a black hole. This will be the essence of the DeDeo-Gott construction. 

Consider a system of two identical conical defects undergone two opposite Lorentz trans¬ 
formations, in the (A 0 , Xl) and (X 3 , X 2 ) planes of the embedding space each: 


^cosli -0 

0 

sinh 0 

0 N 

0 

cosh0 

0 

sinh 0 

sinh 0 

0 

cosh 0 

0 

V 0 

sinh 0 

0 

cosh 0 J 


In the three-dimensional coordinates of the AdS spacetime these Lorentzian transformations 
correspond to SO(2 ,2) isometry transformations. 

It can be shown that in the coordinates of global AdS% these defects move along the same 
circular orbit R = const with a constant velocity, always being at the opposite points of the 
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FIG. 2. (a) A boosted conical defect in AdS 3 . Faces of the wedge are deformed, and the identifi¬ 
cation occurs between points with different time coordinates (in the centre of mass frame). Here 
a = \/ 37 r, ijj = 1. (b) The DeDeo-Gott time machine. 


orbit, Figj^b). A conical defect sits at R = 0 in its rest frame, i.e. 

Xq = cost, Xi = 0, X 2 = 0, X 3 = sinf. 

In the boosted frame 

Xq = cosh 0 cost, X\ = sinh A cos t, X 2 = sinh ?/; sin t, X 3 = cosh -0 sin t, 

hence 

R = \Jxf + X\ = | sinh ^ |, 

X 2 sinh ^ cos t 

cos 0 = — = ——-—— = sgn -0 • cos t. 

R | smh ip | 
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(7) 

( 8 ) 

(9) 

( 10 ) 



























































In the rest frame of a wedge, the points on its edges are identified at equal coordinate 
times. However if we boost it, from the point of view of an external observer this identi¬ 
fication would occur at different times leading to time jumps for a particle moving around 
the conical singularity. In the case of a single conical defect its boost can be regarded as a 
global coordinate transformation of the space-time, which obviously can not cause any new 
physical effects. However, with two defects moving relative to each other, the relative time 
jumps become a physical effect that can not be eliminated by a (proper) choice of coordinate 
system. These time jumps allow for the existence of CTC. 

Existence of closed timelike curves in this space-time can be demonstrated by looking at 
the identification of the edges near the boundary of AdS 3 (at R —» 00 ). We refer the reader 
to [30] for a detailed discussion, here we just quote the result. Speaking in terms of the 
unfolding of the two-conical space-time, when a timclike particle living on the boundary of 
the AdS cylinder hits an edge of one of the two wedges, it undergoes a time and an angle 
jump: 


At = 2 arctan 


Acf) = 2 arctan 


sin(a/2) tanh 7/1 \ 

1 + cos(a/2) tanh ip) ’ 
sin(a/2) 


(11) 

( 12 ) 


tanh -0 + cos(ct/2) 

It can be shown that At + Acft = a. If a > n, the world line of the particle becomes a closed 
timelike curve, and thus the space-time is a time-machine, see Figj3j 

In the regime when this limit is not exceeded, and CTC are not present, the space-time 
has been studied in detail from holographic point of view in [33], but the case of broken 
causality has not been addressed. 

In the next sections we will study geodesic structure of this time machine and explicitly 
show that in presence of the two orbiting conical defects we deal with a highly-nontrivial 
lensing of geodesics, and this reflects on the structure of two point Green’s functions of the 
dual boundary theory. 


III. ENTWINEMENT OF GEODESICS AND CAUSALITY VIOLATIONS 

When conformal dimension A of a boundary operator in AdSjCFT is very high, the 
corresponding two-point Green function can be derived in the geodesic approximation [32]: 

G{A, B) = e ~ ACAB , (13) 
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FIG. 3 . A timelike particle moving along the boundary experiences a jump in time and angle when 
it hits a face of one of the wedges. If At + A(j> > 7r, closed timelike curves become possible. The 
gray strips are the cut out part of the boundary, and the white strips are the physical part of the 
boundary. 

where Cab is the length of a geodesic connecting boundary points A and B. If there are 
more than one geodesic between A and B, they can give additional contributions to the 
propagator. This is the case for the DeDeo-Gott time machine geometry, and here we 
address possible outcome of this in details. 

Consider two arbitrary points A (t 1; <f>i) and B (t 2 , <^ 2 ) located in the physical (unremoved) 
part of the AdS 3 boundary. Having two rotating conical defects in the bulk makes the 
structure of possible geodesics connecting] A and B very nontrivial, so we should find a way 
to calculate their contributions to the two-pont Green’s function G(A, B). Let us shoot a 
geodesic from the boundary point A to the point B. Before it hits the point B it can undergo 
a number of “refractions” on the faces of wedges, winding around either of two conical defects 
clockwise (if it hits the leading face of the wedge head on) or counterclockwise (if it overtakes 
the trailing face of the wedge from behind). For example, schematically a typical geodesic 
1 When A and B are timelike separated we encounter some subtleties caused by the fact that in the AdS 
space-time a timelike geodesic can not reach the boundary. These issues will be commented further on, 
but the general point of view described in this section remains unchanged. 
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may have a structure (see also Fig. [5]) 


A -)• Wf -> -)• W/ -> 5, (14) 

where W^}7 stands for the act of clockwise/counterclockwise winding around the 1-st or the 
2-nd wedge respectively. 

So, formally the Green’s function in the geodesic approximation is given by 

OO 

G(A,B) = E E e -A£(A{Wi-Wn}B) 

n= 0 {Wi-W„} 

where the second sum is taken over all different entwinement structures corresponding to 
the same number of windings, and the first sum is taken over all winding number^} 

It is easy to see that for a given number of windings N the maximal possible number of 
topologically different geodesics is 

n N = 4 ■ 3^- 1 . (16) 


The first winding act can be of four different types. But for each of the next steps, if a 
geodesic wrapped around a conical defect, for example, clockwise, then on the next step it can 
not go in the opposite direction and wrap around the same conical defect counterclockwise. 
It means that in the sequence of windings the winding act W/ can be followed (at least 
hypothetically) by W/, V\?f T , V\?J T , but not by Wf. 

As we will see further, for a given pair of boundary points (A,B) not all sequences of 


entwinements are physically realistic and can contribute to the sum (15). 

To find a proper prescription for the lengths of the non-trivial winding geodesics let us 
discuss in detail a particular example. 

Take a look at Figj4} The length of the presented geodesic is a sum of lengths of its’ 
composing arcs: 


£{A -a WjWtjWJWt, B) = C 


AO\ 


£ 


O2O3 


+ c 


0405 


£ 


0e07 


£ 


OgB- 


(17) 


We will refer to points C) as complementary points, and O* as refraction points. Coordinates 
of the refraction points Oi are to be found from coordinates of A and B , and that can be 
easily done step by step. 

2 This idea of entwinements in holography has been introduced in [[ 34 ] , but there it was related to a concept 
of entanglement entropy “shadows” rather than to subleading contributions to the propagator 
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FIG. 4. Schematic constant time projection of a typical geodesic connecting points A and B in 
the time machine. Red curves are for leading faces of the rotating wedges, and blue curves - for 
trailing faces. The entwinement configuration in this particular case is according 

to the notations introduced in the main text. All shown points in principle can have different 
time coordinates (here we schematically project them down to a single time section, so the curves 
the geodesic is made from should be taken only as an approximate artistic representation). On 
picture (a) the acts of entwinement and identifications are shown explicitly. A and B belong to 
the physical unremoved part of the spacetime, and Oj are the points where the geodesic undergoes 
“refraction” on the wedges. Picture (b) demonstrates the idea of complementary points C t located 
in the removed part of the spacetime. 


Focus on the first refraction on the wedge, Wj■ The points of interest are {A, Oi, C\, C 2 , ()■ 
Here we neglect for a while the second wedge, so we do not consider the point O 3 at all, 
and we treat C 3 as a physical point (not just as a point in the complementary “removed” 
space), see Figg a). These six points can be regarded as a result of boost transformation 
A/ applied to the wedge. We can “unboost” the wedge and find static pre-images of these 


, Ob}- 
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FIG. 5. The left picture demonstrates Wj~ refraction of a geodesic on the first conical defect. The 
right picture demonstrates how it looks like in a reference frame of the conical defect. Again, points 
generically belong to different time slices, and the representation is purely schematic. 

points (see Fig|5](b)). For the “unboosted” points the following relations trivially hold 

C 2 = A - X C 2 = Aj'A - (0, a) = A - (0, a), (18) 

C 3 = Aj 1 C 3 = A fC, - (0, a) = C 1 - (0, a), 

where (0, a) is a boundary identification vector proportional to the angular deficit. Here we 
subtract the identification vector ( 0 ,a) because this particular entwinement is counterclock¬ 
wise. For clockwise W + we should rather add (0,+a). In that case points 0\ and 0 2 are 
intersections of geodesics AC\ and C 2 C 3 with faces of the static wedge. 

Applying the same procedure to the other entwinements, in a generic case we get a system 
of recurrent relations 

C 2 = A 1 (Ar 1 A±(0,«)), (19) 

C 2 j = Aj(Aj 1 C 2 j_ 2 ± (0, a)), 

C 2 N -1 = Ajy(A N 1 B =f ( 0 , a)), 

C 2 j-i = Aj(Aj l C 2 j + 1 =f (0, a)). 
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Here A,- = A/ if the corresponding winding is Wf. A j = A// if the corresponding winding 
is Wfj- In these formulae we pick up the upper sign if W/ /7 , and the lower sign if Wj n . 

Note, that the Lorentz boost we have defined in terms of the embedding space coordinates 
acts non-linearly on the AdSj, points, therefore we can not simply expand the parentheses 


in (19) 


Then for each of the auxiliary arcs C^fcC^fc+i we can derive coordinates of the refraction 
points 0 2 fc, C> 2 k+ 1 , and write down corresponding lengths of the composing arcs. 

Later we will also show that not every formally generated sequence of windings does exist. 

IV. QUASIGEODESICS CONNECTING BOUNDARY POINTS 


To discuss causality properties of the dual boundary QFT, we will in particular need to 
consider boundary points with timelike separation. The conceptual problem we unavoid¬ 
ably encounter here is the absence of timelike geodesics connecting points on the conformal 
boundary of AdS. The equations for timelike geodesics can be derived from the following 
Lagrangian for a massive particle in AdS : 

R 2 


(1 + R 2 )t 2 + 


+ R 2 d) 2 = -1. 


1 + R 2 

Such a particle has two conserved momenta: 

£ = (1 + R 2 )t 1 
J = R 2 j>. 

Substituting them into the Lagrangian we obtain the radial equation of motion: 

£ 2 R 2 J 2 _ 

t -)9 -*-5 


1 + R 2 1 + R 2 R 2 

i? 2 =-(i+« 2 ) (i+fD +?■ 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


Clearly the right hand side of the equation turns negative as R —> oo, and thus no real 
solution to this equation can exist. 

In the case of a stationary spacetime the obstacle could be easily surmounted by per¬ 
forming the analytic continuation of the metric to the Euclidean signature, calculating the 
Green function in terms of Euclidean lengths of the geodesics, and making the inverse Wick 
rotation back to real time. However in our case we deal with a spacetime that is not only 
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non-stationary, but wich has no good global notion of time. Hence we are forced to stick to 
the Lorentzian time. 

The way to implement the geodesic approximation for timelike separated boundary points 
in the single Poincare patch has been suggested in [35]. Let’s turn for a second to the single 
patch of the AdS 3 spacetime, covered by the Poincare coordinates: 

ds 2 = — r 2 dt 2 + + r 2 dx 2 . (25) 

Again, a massive bulk particle has two kinetic invariants: 

E = r 2 t, (26) 

J = r 2 x , (27) 


but now we have two different classes of spacelike geodesics. 


• For J 2 > E 2 : 

r( A) = VJ 2 ~E 2 cosh A 

< x(A) = x 0 + J 2 J _ E 2 tanh A (28) 

t( A) = t 0 + j2 e e2 tanh A 

• For E 2 > J 2 -. 

r{ A) = VE 2 - J 2 sinhA 

< x(A) = x 0 — E 2 J _j 2 coth A (29) 

t( A) = t 0 - e e _ j2 coth A 


We will be interested in the geodesics of the second kind. As A = 0 these geodesics approach 
the point r(0) = 0, which is the Poincare horizon of the half-AdS chart. Regarding the 
horizon as a single infinitely far point (as in the theory of complex functions), we can consider 
two disconnected spacelikc geodesics possessing the same kinetic invariants E and J, but 
emerging from two different timelike separated boundary points A(tA,XA ) and B(tB,Xs), 
as two branches of a single geodesic reaching the spatial infinity and returning back to the 
boundary. The length of such a geodesic will be divergent not only as r —* 00 (the standard 
holographic UV divergence), but also as r —> 0, but this can be cured by an appropriate 
renormalization^} The resulting expression for the renormalized length of the geodesic is 


3 For details see App. B and D of [35] 
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simply 


jC = In ((t B - t A ) 2 - (x A - x B ) 2 ) , (30) 

which gives the correct answer for the two-point correlation function of (1 + l)-dimensional 

cftE 

1 


G(tA,x A ;tB,x B ) = e Ac = 


(( t B ~ t A ) 2 - {x A -x B ) 2 ) 


A ’ 


At > Ax 2 . 


(33) 


In the global AdS space-time the Poincare horizon has no special physical meaning, but 
we can still try to generalize this procedure to this case. 

The boundary held theory now is defined on S 1 x M 1 spacetime, and the two-point Green 
function that we must be able to reproduce via the geodesic approximation has the form 
: 


G(t A , (j) Al is-, 4 >b ) — 


(34) 


i cos(i B - t A ) - cos (4>b - 4>a) I a 
Note that this function is periodic both in angle and time. While the angular periodicity 

is obvious by construction, periodicity in time emerges because of the finite size effects: an 

excitation created at some point in space and time starts dissipating, but later recollects 

and revives due to the spatial periodicity. 

We will need a function that defines angular separation between boundary points while 
properly maintaining the rotational invariance of the system. For instance, given two angular 


coordinates <j >2 = an d <t>i — j , the difference between them along the shorter arc is 


7r 


O(02,0i) = -2 t 


(35) 


Thus we should use 

-D(01, 0 2 )= mod (02 - 01 + 7r, 27 t) - 7T. (36) 

Analogously, for the arithmetic average of two angular coordinates (that provides a point 
exactly at the middle of the shorter arc between 0i and 02 ): 


£(0i, 0 2 ) = 2 


+ 02 — 2nd 


COS ( - (01 - 02) 


4 If Ax 2 > At 2 , the renormalized length is 

C = In (~{t B - t A ) 2 + (x A - x B ) 2 ) 


and the full Green function is 

G(t A , x A \t B ,x B ) = 


1 


\{t B - t A ) 2 - (x A - x b ) 2 \ a ' 


(37) 


(31) 


(32) 
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Note that the Green function (34) has a symmetry: 


G(t A , <j>A] ts, 4>b) — G(t A , tB + X <t>B + 7r) , (38) 

where points on the r.h.s. can be spacelike separated while points on the l.h.s. have timelike 
separation: 


(ye — t A ) 




This symmetry can be used to construct a disjointed spacelike geodesic, with two branches 
reattached at the Poincare horizon, connecting timelike separated points. 

If we represent the global AdS 3 space-time as a cylinder, the Poincare horizon consists of 
two planes cutting the cylinder at 45°. The orientation of the planes (as a rigid construction) 
can be chosen arbitrarily. Then consider a spacelike geodesic emerging from boundary point 
A (see Fig@), and terminating at boundary point B*. Somewhere in the bulk it has a 
turning point P\ where its radial coordinate R* = P(0) is minimal. 

Since we are free to choose the location of the Poincare horizon, we can always orient it 
in such a way that the AB* geodesic intersects it at the turning point Pi. This point Pi can 
be identified with a point P 2 = Pi + ( 7 r t , 71 ^, Or), located on the other cutting plane. The arc 
PiP* can be then rotationally translated to this point: Pi —y P 2 , and then B* —y B, where 
B = B* + (7r,7r). Hereafter we will use disjointed “quasigeodesics” of this APiP 2 B type to 
connect timelike separated points. 

The explicit analytic expression for the quasigeodesics can be derived in the following 


way. Let’s take the Poincare chart geodesics (29), and rewrite them in coordinates of the 


embedding spacetime. The corresponding coordinate transformation is given by: 


X° =^ + 1 + X 2 -* 2 ). 

X 1 = xr, 

x2 = ^- 1+l2 - f2) ’ 
X 3 = rt. 


(40) 

(41) 

(42) 

(43) 


Substituting here (29), and recalling the relations between the global and embedding coor- 
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t = 7T 



FIG. 6. A spacelike quasigeodesic connecting two boundary timelike separated points A and B. A 
spacelike curve emerges from point A and reach the Poincare horizon at point P\. Then it jumps 
to a mirror point P 2 = Pi + (pt- vr^. 0/?) on the second plane of the horizon, and proceeds further 
to B. The arc P 2 B is a rotation of P\B* by 7 r. 


dinates ([3]), we obtain 


y/l + R 2 cos t 
Vl + R 2 sin t 
Rcoscp 
R sin (j) 


[-1 + E 2 - J 2 


J 2 
cosh A, 


2 y/E 2 ^- 
E 

VE 2 - J 2 

J 

y/E 2 - J 2 
(-1 ~E 2 + J 2 
2 y/E 2 - J r 


sinh A, 


cosh A, 


sinh A. 


(44) 

(45) 

(46) 

(47) 


These can be solved to give us the embedding of the Poincare chart spacclike geodesic into 
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the global AdS: 


t( A) = arctan 


2 E 


-1 + E 2 - J 2 


coth A ) + to, 


'_^ _ j ^2 _j_ j 2 

0(A) = arctan ( -—-tanh A ] + 0o, 


2 J 


R( A) = 


J 2 


, 2 x (-1 ~ E 2 + J 2 ) 2 . 2 . 

cosh 2 A + -—^-sirilr A . 


(48) 

(49) 

(50) 


E 2 - J 2 . 4 (E 2 - J 2 ) 

This quasigeodesic already has a 7r-jump in time at the turning point A = 0, and as explained 
before we also need to adjust the discontinuity in angle: 

_1 _ j^2 _i_ j2 \ 

-—-tanh A J + tt9( A) + (j) o . (51) 

From now on we will omit the tilde. 

The kinetic invariants can then be expressed in terms of the boundary coordinates: 


E = 


J = 


Sill 


t 2 ~tl 


sin -cos^’ 


cos 


2 2 
D(<f>2—Tr,4>i) 


sin D(02 7’ 0l) - cos ^ 


The integration constants t 0 and 0o bi (49) can be represented as: 


— 2 (^i + ^2) , 0o — S(0i, 02 — tt) • 


(52) 

(53) 

(54) 


Inverting equation (50) we obtain dependence of the affine parameter on the radial coor¬ 
dinate: 


A (R) = ± arcsinh 


4 (E 2 - J 2 )R 2 - AJ 2 


(55) 


(-1 - E 2 + J 2 ) 2 + 4J 2 ’ 
where the minus sign is taken on the first branch of the geodesic or quasigeodesic (i.e. before 
the turning point, - when particle moves away from the boundary), and the plus sign is taken 
on the second branch (when particle moves towards the boundary). This function can be 
used to define the geodesic length, which is simply 


C(R u R2)= \±(R 2 ) - \±(Ri) , 

for two points with radial coordinates R\ and i?o- 

Note that for a geodesic connecting two boundary points the length is divergent: 


(56) 


jC = lim (A +(R) — A_(i?)) = 2 lim 


R-* 00 


A(E 2 - J 2 )R 2 - AJ 2 

o y (-1 - e 2 + j 2 ) 2 + aj 2 


= 00 


(57) 
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and needs to be renormalized. The natural way to do it is to subtract the parameter 
independent divergent part, and define the geodesic length as 

16 (E 2 - J 2 ) 


C ren = lim (A_| .(R) — X-(R) — 2Ini?) = In 

R—>oo 


(58) 


E* - 2E 2 (—1 + J 2 ) + (1 + J 2 ) 2 , 

The argument of the logarithm is always positive for quasigeodesics connecting timelike 
separated points, but can be less than 1. It means that the geodesic length after the 


renormalziation in principle can be negative. Using (58) along with (52) we obtain the 
correct result for the retarded Green function: 


G c (t i, 0 2 ) — e 


— AC r 


(cos (t 2 - ti) - COs(0 2 - 0l)) A 
where (t 2 - U) 2 > D(0i, 0 2 ) 2 , 


(59) 


where Z)(0 1 , 02) is the function introduced in (36). The possible negativity of the renormal¬ 
ized geodesic length is the reason why singularities of the correlator can be captured in the 
geodesic approximation. 


Here we must pause for a second and stress Lorentz non-invariance of (59). We define 


50(2,2) isometries of AdS% in terms of the Lornetz boosts of the embedding id 2,2 } space 
(J6|. If we take two boundary points A and B, and act on them with a bulk isometry 
transformation A of this type, we will observe that it does not preserve the Green’s function 


(59): 


G C (AA,AB)^G C (A,B). 


(60) 


It is the fundamental difference between holography of a Poincare chart and holography of 
global AdS. In the first case the bulk isometries induce Lorentzian boosts on the boundary, 


so the Green function of a dual boundary field theory is a relativistic invariant object (33). 


In the second case the isometries rather act as conformal transformations leaving the Green 
function covariant, i.e. invariant up to some coordinate dependent scaling prefactors. 

In the holographic language this is encoded in the fact that the renormalized lengths 
connecting boundary points are dependent on the choice of the reference frame. Below when 
we consider the DeDeo-Gott time machine geometry, we should be especially careful about 
this, since the geodesics there are combinations of Lorentz invariant and non-invariant terms 
as, for example, in ( Jl7| ). The proper way to deal with it is explained in the first subsection 
of Sec. El 
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In the next section we will analyze lensing of the quasigeodesics on the conical defects 
and calculate the Green function of the dual held theory in presence of the closed timelike 
curves in the bulk. 

V. THE TWO POINT GREEN’S FUNCTION 

A. The algorithm 

In Sec. |m]we have discussed the general idea of using the geodesic approximation to com¬ 
pute the boundary Green’s function for the DeDeo-Gott geometry. Now we will formulate 
an exact algorithm for that. 

• Introduce coordinate system on the unfolding of the double-cone space in such a way, 
that the physical (unremoved) part of the boundary consists of two stripes covered by 
coordinate intervals: 


t G (— 00 , 00 ) , 


(61) 



• Fix two boundary points A and B. For simplicity we can choose A = (0, — §). 

• Fix the total number of windings N that a geodesic of interest undergoes on the way 
from A to B. In our simulations we will not go beyond N — 4, because the higher-order 
contributions to the Green’s function are highly suppressed. 


• For the given A, B , and the number N, generate all possible 4 ■ 3^ 1 sets of the 


complementary points {Ci, • • • C 2 n} corresponding to different sequences of windings 
{W 1 , ■ ■ ■ Wm}. The (quasi)geodesics then consist of N + 1 arcs AC\, C 2 C 3: ..., C 2 nB, 
each of which is just a (quasi)geodesic curve in empty AdS 3 . 

• Impose that each of the “odd” complementary points C 2l +\ belongs to the causal future 
of the previous “even point”: C^+i >- C 2 i. 

Let us elaborate on what the reason to do so is. The most clear question we can 
ask is whether causal propagation of a signal from the future to the past is possible. 
To define the dual retarded Green function in presence of the CTC in the bulk, we 
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should recall that evolution of a particle moving in the bulk of AdS can be split 
in two parts: “physical” continuous motion along a timelike or a spacclike geodesic, 
and “topological” time jumps caused by entwinement around the conical defects. In 
the holographic language geodesic branches C^C^+i correspond to the continuous 
evolution, and C^+i —> € 21+2 identifications - to the time jumps. In absence of the 
closed timelike curves a signal could causally propagate from A to B if B belongs to 
the future light cone of A: B >- A. A natural generalization of this prescription for the 
time machine case is to impose that this should hold true for all “physical” segments, 
he. C* 2 j+i C* 2 j V i. 

• For each of the causal quasigeodesics, solve for the intersection points {Oi, • • ■ O 2 N}■ 
The easiest way to do this is to transform for each winding back to the rest frame of 
the corresponding wedge. For example, if branch C^C^j+i intersects first the trailing 
face of the 2nd wedge, and then the leading face of the 1st wedge, we perform a Lorentz 
transformation of the branch to the 2nd rest frame, then untwist the wedge by a simple 
coordinate transformation <j>' = (f> — t, such that angular location of the face remains 
still in these co-rotating coordinates, and solve the equation 

($Lc 2i+1 ( A 2 i) ~ tc 2 iC 2 i+ 1 ( X 2 i)) mod 2 n = (j>' Tn mod 2vr . (62) 

Then we repeat the procedure in the 1-st rest frame: 

(<^c 2l c 2i+1 (-W) - tc 2 i c 2i+1 (A 2 i+i)) mod 2vr = <f>' Li mod 2vr . (63) 

• Make sure that all these equations have real solutions (otherwise discard the geodesic). 

• Make sure that if a branch is not expected to intersect other faces within the physical 
region of the space, it actually does not (fake intersections within the removed part of 
the unfolding are allowed). In other words, if an arc OjiOji emerges from the face Lj 
and terminates at the face Tji, it should not have intersections with L n and Tj. 

• Calculate the lengths of all inner segments of the geodesic (C> 2 i 02 i+i). They are finite 
by construction and equal to 


L 


02iC>2i+l 


'2i+l 


— A 


2 i 


(64) 
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FIG. 7. A quasigeodesic with winding number N = 2 connecting boundary points A and B. 


• Renormalize the lengths of the boundary segments AO\ and C> 2 nB as they are diver¬ 
gent: 


£.40, = + *i. (M) 

Lo 2n b = 2 l c?nB - -W - (66) 

• Calculate the renormalized lengths of Lao x and Lo 2N b in the original frame. As 
mentioned in the previous section, the renormalized lengths are not Lorentz-invariant. 
So, while we are free to constantly switch between different reference frames in order to 
calculate lengths of the finite inner segments C^C^+i, the renormalized lengths of the 
two boundary segments must be calculated in the original frame where we define the 
Green’s function. In our case it is the “centre-of-mass frame”, where the two conical 
defects are symmetrically boosted. 

• Finally calculate contribution of the geodesics to the Green function: 

G(A,B)=y e -“', (67) 

k 

where the index k runs over the set of geodesics that satisfy aforementioned conditions. 
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In accordance with the described algorithm we subsequently account for contributions 
from higher winding numbers starting with N — 2^J In other words, we formulate a kind of 
“perturbation theory” with the number of entwinements as a control parameter. 

One property of this series expansion must be comment on. Each geodesic contributes 
to the Green function exponentially: 

e" Ai ™. (68) 

For higher windings the number of internal segments C^jCWi the geodesic grows linearly 
in N, and so does its renormalized length L ren . Therefore the corresponding contribution to 
the Green function is exponentially smalj^J On the other hand the total number of possible 
topologically different geodesics scales as 

n ~ 3*" 1 , (69) 

i.e. grows exponentially. Therefore in principle these two effects can compete and we can 
not say a priori that the higher order contributions to the Green function are suppressed, 
and the sum over entwinements is convergent. If not, this could mean that our setup is 
unstable and undergoes a Hagedorn like transition. 

However there are three different reasons for it not to happen. Firstly, by no means all of 
these 4 • 3 N ^ 1 winding configurations satisfy the causality condition: C^+i >- CV Secondly, 
even if the causality condition for the set of complementary points is satisfied, the geometric 
structure of the geodesics becomes more and more complicated as the number of windings 
increases, and it becomes hard to force a geodesic curve to undergo the concrete sequence of 
windings (it is easy to see on Fig|8](b)) . Finally, the “decaying” exponent has a conformal 
dimension as a knob, so at large enough A it dominates over the “growing” exponent. 

Another way to understand convergence of the series expansion in all orders relies on a 
simple and general argument. Consider a germ of all possible quasigeodesics emerging from 
point A. The Erst segment of a generic quasigeodesic curve hits the boundary at some point 
Ci, first in the sequence of complementary points {Cj}. If we go along the curve further, 
we will obviously see that it is defined uniquely up to the final moment when it reaches the 
physical part of the boundary at point B. Thus, for a fixed initial point A, for each of the 

5 For negative times N = 1 windings do not contribute as they are due to the lensing on a single conical 

defect, that obviously can not lead to time travelling. But for positive times we take them into account. 

6 In a generic case when L ren > 0. 
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FIG. 8. (a) An example of a physical N = 4 geodesic contributing to the Green function, (b) An 
example of physically impossible winding: for a given sequence of entwinements a geodesic can not 
be fit in the unremoved part of the spacetime. 

“first-in-the-sequence” complementary points C\ the final point B is defined unambiguously. 
Now, as we have already emphasized, the singular contributions to the Green function come 
at the points where the renormalized geodesic length is infinitely negative, C ren = — oo. It is 
possible if and only if C\ is located exactly on a generatrix of the light cone emerging from 
point A, or C 2 N is located on the generatrix of the light cone of point B. Quasigeodesics 
having the complementary points C\ and C 2 n right on the corresponding light cones form 
a zero measure subset among all possible quasigeodesics. Thus the set of boundary points 
where the Green function is infinite is also a zero measure subset of the boundary spacetime. 
Everywhere else the Green function is finite and well-defined. 

B. Phenomenology and discussion 

We are now ready to implement our computational algorithm for the DeDeo-Gott geom¬ 
etry. As shown in Sec. [TT] closed timelike curves in the spacetime are present when the total 
angular deficit is more than 2n. For concreteness we impose oljjt = \Z3n, and the boost 
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FIG. 9. The Green function is sourced at point A with coordinates (0, —tt/2). Blue lines depict the 
light cone emerging from this point. Any timelike line would cross the removed region and enter 
another strip (red dashed line). In order to avoid formulating the boundary field theory on both 
strips simultaneously we calculate the Green function on a timelike line very close to the light cone 
generatrix. 


rapidities 'i/jjjj = ±1. In their corresponding rest frames (in the co-rotating coordinates) 
the locations of the edges are taken to be 


= «/2 0ri = 

4>'l 2 = a/2+ 7i 4>' T2 = -a/2+ 71. 


(70) 


It is more convenient to calculate the Green function also in the co-rotating coordinates: 


G° r (tl, t2, </4) — G(t\, 01 — t\] t2, 02 — G) ■ 


(71) 


For simplicity we will mostly study the Green function on a one-dimensional timclike line 
passing through the point A. We should be careful here. Any timelike line originated in 
one physical strip crosses the cut out region and enters the second strip. To formulate a 
quantum field theory on both strips simultaneously is possible yet tricky due to the fact 
that on the unification of two parts of the boundary time can not be globally defined. To 
avoid this difficulty we will consider the Green function on a timclike line in a close vicinity 
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of the generatrix of the light cone. Then in a large range of times we will stay within one 
strip of the boundary. 

In other words, the object we will attempt to evaluate is (in the co-rotating frame) 

G cr (0, —vr/2; t, — 7t/2 + et) , e<l. (72) 

We have performed the numerical calculation of the retarded Green function for negative 
times t < 0 up to N — 4 order, and for positive times t > 0 up to N — 2. 

Let’s discuss firstly the analytic behaviour of the Green function at negative times, - how 
the quantum particle behaves travelling back in time. A naive expectation would be to think 
that the Green function decays as t —>■ — oo, and it is partially true. However as we can see 
at not very large negative times the function develops a number of non-trivial features, - 
peaks which we can interpret as the “most probable” regions of times the particle can reach 
using the time machine. The corresponding results are present on FigflOj As an illustration 
we also provide a two dimensional plot for the leading N = 2 winding at negative times, 

Fi® 

The origin of these peaks can be traced back to the fact that renormalized length of 
a geodesic can be negative. Generically at small conformal weights N = 2, N = 3 and 
N = 4 contributions are commensurate, but already at A > 2, higher entwinement terms 
are getting suppressed as compared to iV = 2. However at specific points, where C ren < 0, 
the corresponding contributions to the Green function are getting enhanced in the large A 
limitQ forming a sharp peak. For instance, N = 2 set of geodesics contains such a curve 
around £2 = —1-9, and N = 4 set has a special point at f 2 — —1.45, FigjTO} 

We have not performed numerical simulations for N > 4, but we can not exclude that 
such negative length curves can appear also at large N. The geodesic length is defined by 
lengths of internal segments (always positive) and lengths of the two boundary segments 
(that in principle can be negative): 

N 

£ ren r*ren . nren . \ ^ r 

AB ^AOl ' *~'C>2 nB ' / ; / ~'02i-202i-l ■ V'H 

1=1 

N 

!f A 01 + £o 2 N b < °> and Wao, + £oZb\ > s £o 2i _ 2 o 2i _ i; the Green function will get 

i=l 

a contribution that does not vanish in the large A limit. For a large number of internal 
7 Strictly speaking, the geodesic approximation is reliable only in this limit. 
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G(0, — 7r/2; t, (1 — e)t — 7t/2) 



FIG. 10. (a) N = 2, (b) N = 3 and (c) N = 4 contributions to the retarded Green function at 
negative times at A = 1.5. Discontinuities of the curves are artifacts of the geodesic approximation, 
(d) The retarded Green function at A = 1.5 (N = 2, N = 3 and N = 4 contributions are added up). 
For the large conformal dimensions peaks are enhanced, not suppresses, and we can see revival of 
the particle at moments preceding the excitation of the Green function. A not very large conformal 
dimension is chosen for convenience of presentation. Here e = 0.1. 

segments it is not likely, but neither is impossible: while all internal lengths are finite, the 
renormalized negative lengths might be of an arbitrarily huge absolute value: 

L r lo,< 0, I^moJ » !. (74) 


thus dominating over positive contributions. 





























FIG. 11. N = 2 contribution to the retarded Green function at A = 1.5 at negative times in two 
dimensions. The Green function is sourced at the red point (0, — 7t/2). We construct the theory 
only within one of the two strips of the boundary. Sudden break of the function signalizes that 
some regions of the spacetime are unattainable for the N = 2 quasigeodesics. 

In the case of a large conformal dimension it would mean that, if we were able to sum 
up contributions in all winding orders, the resulting Green function would have a shape of a 
comb with a number of peaks (in our calculations we discovered two of them). These peaks 
play a role of “pit stops” for a particle travelling in time, - they form a set of easily reachable 
coordinates in time. Hence we deal with specific “negative time” revivals. 

At positive times we have also discovered interesting features of the Green function. In 
the case of plain AdS 3 geometry the dual light-like Green function (shifted away from the 
singularity) is decaying in time. I 11 presence of the causality violating conical defects we 
detected a new peak of a high weight, Figjl2| signaling a revival of the excitation. 


VI. PHASES OF THE BOUNDARY FIELD THEORY 

We have calculated the Green function numerically up to N = 4 entwinements for the 
time machine geometry with a = \Z3n and i/j — 1. However it would be interesting to study 
how the properties of the Green function change upon changing the strength and rapidities 
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G( 0, — 7r/2; t. (1 — e)t — 7r/2) 



FIG. 12. The retarded Green function at conformal weight A = 1 at positive times. The plot 
demonstrates how the non-trivial N = 2 windings modify the original Green function ( N = 1 
does not contribute when we consider the boundary theory within a single strip). The offset plot 
represents the N = 2 winding contribution separately. We have made the Green’s function timelike 
to avoid dealing with the light-cone singularity. Here e = 0.1. 

of the conical defects. We constructed the leading order N = 2 contribution to the retarded 
Green function at negative times for a G (1.17T, ...1.957t), and 0 £ (0.1, ...1.5) with stepping 
A a = 0.057T, A0 = 0.05., paying special attention to the location in time and and strength 
of the revival peak. 

The results can be schematically summarized in a form of a phase diagram, Figjl3j 

• If for a given value of a the rapidity 0 is not large enough to prevent the system from 
collapsing; the DeDeo-Gott geometry is forbidden (blue). 

• If for a given value of a the rapidity 0 allows for the existence of the DeDeo-Gott time 
machine, but still not very large, we clearly see the effect of revival, and the peak is 
sharper the closer 0 is to the lower bound (yellow). 

• If the rapidity is too large, the causality is violated, but excitations just decay and do 
not revive at negative times anymore (green). 

• At very small values of a the retarded Green function does not exhibit any non-trivial 
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1.17T 1.37T 1.57T 1.77T 1.97T 

a 

FIG. 13. The four different “phases” of the boundary field theory. Blue spots depict the region 
of forbidden geometries, yellow spots are for the phase of negative time revivals, green spots 
form the region where the retarded Green function is non-zero at negative times, but does not 
exhibit reviving peaks in the large A limit. Red spots are where at the leading N = 2 order 
the boundary field theory retarded Green function does not demonstrate causality violation (i.e. 
G cr ( 0, —7 t/2 ; t, — 7 t/2 ) = 0, t < 0) despite the presence of the CTC in the bulk. Everything is based 
on the numerical simulations of the leading N = 2 contribution to the retarded Green function. 
We expect higher order corrections to change the diagram qualitatively, but not quantitavely. 

features at negative time even in presence of the closed timelike curves (red). However, 
this feature is likely just an artifact of N = 2 approximation, and we do not expect it 
to be there for higher windings. 

The profiles of the Green’s function at negative times are presented on Figjl4)(a,b) for 
a = 1.57T and a = 1.77T respectively. The fact that revivals are seen only at not very large 
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G(0, —tt/ 2; i, (1 - e)t - tt/2) G (°: t, (1 - e)t - tt/2) 




FIG. 14. The negative time N = 2 contribution to the retarded Green’s function at (a) a = 1.5, 
ip = 0.45, 0.50, 0.55, 0.60, and (b) a = 1.7, ip = 0.75, 0.80, 0.85, 0.90. Here A = 1. On each of the 
plots the first three peaks are getting stronger and sharper as A —» oo (the yellow region on the 
phase diagram), while the fourth one is getting suppressed in the same limit (the green region), 
so we do not consider it to be an actual revival of a non-causal excitation. Counterintuitively, the 
weight of the Green function drastically decreases when we increase a or/and ip. 


rapidities (and the effect is stronger as closer ip to its minimal possible value) is surprising 
and contrasts to how causality is broken in the bulk. The stucture of CTC is defined by a 
and ip , and the time jumps become stronger as the angle defects and rapidities are increased. 
Thus we rather should expect that for high ip the time travelling along the CTC is more 
efficient in the sense that amplitudes of the classical free Green’s function defined on the 
boundary are getting enhanced as a ip grows. In the interacting holographic dual field theory 
the retarded Green’s function is damped for larger a and ip , so we can claim that causality 
in the boundary field theory is broken mildly as compared to the bulk. 

Another interesting feature of this system is that while the overall weight of the Green 
function drastically decreases when the rapidity ip is taken away from the “forbidden region” 
on the diagram, the actual past time penetration depth (i.e. the deepest reachable point at 
negative times where G cr (0, — 7t/2; t, — 7t/ 2) yf 0) increases (though very moderately), and 
this is in agreement with the “naive” intuition. 
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VII. CONCLUSIONS 


In this paper we have analyzed properties of a two-point Green function in a (1+1)- 
dimensional held theory dual to the DeDeo-Gott time machine geometry. Using the geodesic 
approximation we have shown that AdS/CFT is capable of describing a quantum held theory 
when causality is violated, and have shown that the corresponding boundary propagator has 
remarkable features. We discovered that in presence of closed timelike curves in the AdS 
bulk a causal propagation of an excitation from the future to the past is possible on the 
boundary, and the retarded Green function exhibits peaks at certain negative times. At 
positive times analytic structure of the Green function also changes, and new singularities 
arise. 

Surprisingly, we have found that as we increase the strength of the conical defects a and 
the rapidity % i, the causality violation in the dual held theory is getting milder in the sense 
that the weight of the retarded Green’s function at negative times decreases. 

Contra to the previous results on the dynamics of physical systems in time machine 
backgrounds [23H27j our calculations have demonstrated that sometimes it is possible to 
define evolution of an interacting theory in a time machine without imposing any additional 
self-consistency constraints. Despite the explicit non-causality the Green function does not 
have any uncontrollable pathologies. 

Our considerations leave a number of open questions. First of all, we have to understand 
how to interpret the boundary state dual to the DeDeo-Gott geometry, - whether this quan¬ 
tum state is pathological or just exotic yet physical state. From the boundary point of view 
a single conical defect, if its angular deficit is a = 27t( 1 — 1/IV), can be thought of as state 
created by a non-local twist operator in a conformal field theory [33]. But what it means to 
have such an interplay of two independently boosted defects has to be clarified. 

Another thing we have not touched on in the paper is the entanglement structure of the 
boundary state. We focused on the properties of the retarded Green function, and thus 
analyzed the timelike quasigeodesics. However, even below the a = n threshold, when the 
CTC are not present in the system, due to the lensing it is possible to connect timelike 
separated boundary points just by standard continuous spacelikc geodesics. If a certain 
generalization of the Ryu-Takayanagi conjecture 133 is true in this case, it would mean 
the boundary state is timelike entangled [38]. Possible physical outcomes of this fact is an 
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interesting direction for future research. 
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